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Abstract 

We determine the damping rates of slow-moving photons in next-to-leading order hard-thermal- 
loop perturbation of massless QED. We find both longitudinal and transverse rates finite, positive, 
and equal at zero momentum. Various divergences, light-cone and at specific momenta, but not 
infrared, appear and cancel systematically. 
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Hard thermal loops (HTL) have been used recently in different contexts. For exam- 
ple, to discuss CP asymmetries and high-temperature leptogenesis [fj, to determine the 
three-loop HTL contributions to the thermodynamic functions of a quark-gluon plasma ^ 
and in finite-temperature QED [3], to calculate the electric and magnetic properties of the 
quark-gluon plasma |4|, to find the quark number susceptibility at finite temperature and 
zero chemical potential More theoretically-oriented work covers trying to improve the 
HTL perturbation expansion 6[, finding a recursion relation for special one- loop diagrams 
that translates into an equation for their generating functional [7| , calculating hard thermal 



loops for a spatial 't Hooft loop in the deconfined phase of a gauge theory 



8j, studying 



the high-temperature behavior of n-point thermal loops in static Yang-Mills and gravita- 



3, 



lookin 



g at a 



11] , and 



tional fields [9], studying the real-time formulation of hard thermal loops 
HTL-resummed improved ladder Dyson- Schwinger equation and gauge invariance 
extending the formulation to M = 4 supersymmetric Yang-Mills theories jl^ . 

Hard thermal loops have been the answer to early problems encountered in standard 

n n 

loop-expansion in high-temperature QCD [131 ] . In a line of works [14J, we have looked at the 
behavior in the infrared of HTL-fully-dressed one-loop dispersion relations of slow-moving 
longitudinal and transverse gluons, quarks, fermions in QED, and quasiparticles in scalar 
QED. The present work finishes this line by looking at the one-loop HTL-fully-dressed 
damping rates for slow-moving longitudinal and transverse photons in small-coupling finite- 
temperature massless QED. These are found finite, positive, and equal at zero-momentum. 
Some of the issues encountered in the context of finite-temperature QED can be found in 
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In the imaginary-time formalism, the quasiparticle momentum is P^ = (po> p) with P 2 = 
Pq+P 2 . The Matsubara frequency po is equal to 2mrT ((2n + 1) ttT ) for photons (fermions), 
n an integer. Real-time amplitudes are obtained via the analytic continuation p = —iu^ + 
+ , with the real energy of the quasiparticle. We assume a regime of small coupling 
constant e so that the temperature T is the hard scale, eT the soft scale. The poles of 
the HTL-dressed photon propagators A^(— iw,p) determine the dispersion laws oj^tip) at 
lowest order in the coupling (I for longitudinal, t for transverse). We use the Feynman 
gauge. The lowest-order dispersion relations write: 

-uj +p 2 + 8Tl e (-iuj £ ,p) = 0; -oo 2 +p 2 + 5U t (-ioo t ,p) = 0. (1) 



2 



The quantities 5He and 5H t are hard thermal loops, given by: 



SU,(-iu;,p) = ^-1 [^-Qo(( £ ) 

p 6 \0J Z J lu \oj/ 



5U t (-iu,p) = 3wq 



p 

00 



(2) 



with Q (x) = | log ( J and u = eT/3 the frequency of the plasma. For soft momenta, 



the dispersion laws, real at lowest order, can be determined analytically: 



Ue (p) — Wn [ 1 H T? p 4 



(p) - ' 1 • tp' 2 Jr/' ; • 



(3) 



with p = p/uiQ. 

The HTL-one-loop dressed propagator D^ v (P) is obtained by adding to the lowest-order 
A~J (P) the HTL-dressed one-loop-order photon self-energy n^(P): 



D-, 1 (P) = A" 1 (P) + IV (P) • 



(4) 



Using the transversality identity Pfjl^ = 0, we can decompose the dressed self-energy into 
longitudinal and transverse components, namely n Ml/ (P) = 11^ (P)P^,+ II t (P)P*„, with Pj^ 
the corresponding projectors, and: 



u t (P) = ^ n 44 (P) ; n t (P) = I (n,, (p) - (p)) . 



(5) 



2 

Sum over repeated indices is assumed and p = p/p. The damping rates j it (p) are the 
negatives of the imaginary parts of the corresponding complex energies, poles of D^ t (P). 
In HTL-dressed perturbation, they are obtained via the relation: 

ImlL^ (-iu,p) 



lt,t (?) 



— 2uj + dSIl( jt (—iu,p) /du 



(6) 



u>=ue,t(p) 



The denominators in this expression are readily calculable, and we have: 



UJ 



1 



7^) = ^[l + 0(p 2 )]lmIT 44 ; 7t (p) = ; 



:P 2 + ^ (P 4 ) 



Im(IIii -piPj-IIij) . 

(7) 

The HTL-dressed self energy n is the sum of two one-loop diagrams and n( 2 ) given 
by the relations: 



11$ (P) = e 2 Tr soft Tr 7 [I^tf - P, -X; P)Af(K)T u (K, -K + P; -P)A f (K - P)] ; 
ng) (P) = e 2 Tr soft Tr 7 [I^tf, —if; P, -P)A / (if )] . (8) 



3 



In these relations, Tr so f t = J d 3 k/ (2tc) 3 and 'soft' means only soft values of the loop 

n 

momentum k are to be included. Also, Tr 7 means tracing over Dirac gamma matrices. The 
amplitude is the equivalent in HTL perturbation of the classic fermionic loop whereas 
is a pure HTL effect. The propagator in ([8]) is fermionic, writing as: 

dVL s $ 



Af (P) = A o ; (P) - SE f (P) ; 52f (P) = m)J 



4tt PS' 



(9) 



The quantity Ao/ is the undressed fermionic propagator and SH a hard thermal loop. The 
quantity rrif — a/1/8 eT is the thermal fermion mass and 5* = (z, s). We can decompose A/ 
along helicity eigenstates: 



A f (P) = - [ 7+p A + (P) + 7 _ p A_ (P)] . 



(10) 



Here, ~f± p = (j ± 27. p) /2 are the two helicity-state projectors and 7^ the euclidean Dirac 
matrices. The two propagators A± are equal to (D =F P s ) _1 respectively, with: 



D (P) = ip --LQ MS ; Ds (P)= v + ^ 
P V P J P 

The HTL-dressed vertices F involved in (j8] ) are as follows: 
T,(K-P,-K; P) = l, + m}J 
F^(K,-K;P,-P) -- 



P VP 



(11) 



4vr PSQS ' 



-2ml 



(12) 



J 7 47r PS* (P + K) S (P — K) S' 

In view of ((Tj), we only need the quantities n 44 , and pipjUij. We write these for the 
first amplitude in (JHJ) as follows: 

e 2 r 

nff = - jTr soft ^ [l + erjkq + 2I + 2el 3 + 2tjl, + 2e V kq/ 

£,77 

+ 2e/ / 3 + 2 V I I Q + 2e V I 3 I q + J 2 + II + J 3 2 + erjkq (/ 2 - J 2 - if) 



A e A„; 



(13) 



H« 



g2 r 

-Tr soft ^ [3 - ET] kq - 2 £ / 3 - 2I + 2£7/kqi - 2t?/ 9 - 4er]I 3q 



+ 2 £ i 2 i 23 + 2ei 3 i 33 + J 2 + i 3 2 + 4 + / 2 2 + 2i 2 3 + I 2 33 + 2r]I 2 I 2q + 2I 3 I, 



+ 2er]I 23 I 2q + 2er]I 33 L 



3q 



'- 22 



T 2 _|_ 2I 2 



23 



I 2 ) 



3q 



(14) 
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e 2 r 
PiPjU-if = — Tr soft ^ |l - er/kq + 2erjxqp 



e.t) 



- 2r/qp/ p - 2exl p - 2I PP - 2erjqp I p3 + 2er]kq I pp - 2er]xl qp 

+ I p + 2eI p I p3 + 2 v I p I qp + 2e V I p3 I qp + I 2 pi - e V kq (/J - J 2 )] A e A^. (15) 

In these expressions, A e = A e (K) and A^ = A,, (Q) with Q = K — P. Also, the 
summation indices (e, rf) belong to the set {(+, +) , (— , — ) , (+, — ) , (0, +) , (0, — )} with 
A (K) = Q (ik /k). Furthermore, x = pk. The quantities lap... are solid-angle integrals 
given in (fTT]) below. For the second amplitude in ([8]), we have similar expressions: 



ni? (P) = - Ili 2) (P) = -e 2 Tr soft ^ [(A; - A e ) J + eA' £ J q - sA £ J 3 ] ; 
frUf) (P) = e 2 Tr soft ^ [(A^ - A £ ) J pp + eA £ J ppq - eA e J pp3 j . (16) 



The quantities J Q p... stand also for solid-angle integrals, defined here: 

2 fdtt s asf3s... 2 [ dn s as frs . . . 

u - = m fJ 17^s^ ; J «*~ = m fJ 17 Ts^ks- (17) 

For both integrals, the indices a, ft, . . . correspond to the directions x, y, z, p, k and q. The 
subscript means 'one' in the numerator of both integrands. 

To proceed, it is best to work out with some detail one typical term. Let us therefore 
consider the following contribution to longitudinal Ilj^, namely the {H — } term in the sum 
$13]). We write: 

n+-(P) = n 5 (P) + n c (P) ; 

e 2 

n s (P) = - jTr soft [1 - z + 2 (1 - z) Jo + 2/ 3 - 2I q ] A + A_; 

n c (P) = -^Tr soft [2/0J3 - 2I I q - 2/3/, + (1 - z) J 2 + (l + z) (J 2 + J 3 2 )] A + A'_ . (18) 

In these expressions, we have replaced kq by z = (k — px) jq. Also, from now on, we take 
Wo = 1 to ease intermediary notation. It is easy to show that the integral I in ( [T7|) is the 
difference b — a where: 

.,, . dQ s as /3s . . . 2 f dfl s as /3s . . . 

<*- =m fJ 17 -psw^> h <*~ =m n i^^sQ^' (19) 

so that, with proper replacements, we get the result: 

e 2 

U S (P) = --Tr soft [1 - z - 4 (1 - z) a - 4a 3 + 4a g ] A + A_, (20) 



and the results: 



nc(P) 
n-(p) 

I^(P) 



u a c a (P) + n™(P) + n^?(P); 

e 2 

-— Tr soft [2a a 3 - 2a a g - 2a 3 a 9 + (1 - z) a 2 , + (1 + 2) (a 2 + a 2 )] A + A_ 

e 2 

--Tr soft [26 &3 - 2b b q - 2b 3 b q + (1 - z) b 2 + (1 + z) (6 2 + 6 2 )] A + A_; 

e 2 

yTr soft [a 6 3 + a 3 &o - a & g - 0,60 - a 3 & g - a 9 6 3 
+a & (1 - z) + (ai&j + a 3 6 3 ) (1 + z)] A+a'_ . 



(21) 



The angular integrals involved in a a can be obtained in a power series in p. In the longitudinal 
case and up to order p 2 , we have the following results: 

9 



a 
ax 



8k 
9 

8k 



a 2 = 0; 
9 

8k 
r 



«3 



-xp + T - (l - 3x 2 ) p 2 + ^1 + rxp + (l5rV - 5x 2 - 5r 2 + 2) p 2 ^j Q + . . 
~P ~ ^3r 2 - xp 2 + + rxp 2 ^ (r 2 - l) Q + • • 



-L / . _ O Q O \ O r O 

1 + rxp H (45r x — lor + 1) p — [r + r xp 

30 



— (l5rV - 5x 2 - 5r 2 + 2) p 2 }Q + ... 



a 



x — tzttT (l — 3x 2 ) p H — — (l5r 2 — 4) (hx 2 — 3) xp 2 



'" 8fc 16fc 
9 

8^ 



80fc 

rx + — (3r 2 x 2 — x 2 — r 2 + l) p + — (5r 2 — 3) r (5x 2 — 3) xp 2 



Qo 



k p 

Ouq = — (2 3 — — dp. 

q q 



(22) 



In these expressions, r = ik /k and Q = Qo( r )- The angular integrals 6 a in ( !T9~]) can 
be obtained from these with suitable replacements and further expansion. Plugging these 
expressions in fl20l) . we obtain 1 : 
q _ ^ 



n s (P) = e 2 p 2 Tr soft - 



16A; 3 



[(9 - 9kr - 2k) - (8k 2 r - 8k 2 - 9) (k + kr - 1)] A + A_. (23) 



In this expression, terms involving the product of three functions that will necessitate later 
a spectral decomposition, heavy on the extraction of the imaginary part, have been worked 



1 In these longitudinal calculations, the zeroth-order contributions cancel as they should and hence are not 
shown. Only O (p 2 ) contributions are kept and displayed. In view of this will yield 7^(0). 
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out using the following identity: 

QoA+ = 



1 



8k" 1 



8k 



(24) 



r-1 9 J T ' 9(r- 1) 

Similar work is done for lie. The aa and 66 contributions, making suitable changes in the 
latter, yield together the result: 

(1 -x 2 ) 



Tl a c a (P) + II™ (P) = eV Tr 



256A; 4 



x 



[(-81 (At - l) 2 + 18 (8A; 2 - 8k 2 r + 9) 
1 - kr - k) (1 - kr) - (8k 2 r - 8A; 2 - 9) 2 (k + kr — l) 2 ) A + a'_ 
. +72k (k + kr- if (r - 1) Q A_ 



(25) 



However, the third contribution requires more work now that a and h integrals are side 
by side. Still, this can be done and we find: 

;i-x 2 ) 



W$ = eVTr soft V 256fc2 ; (8k - 8kr - 8k 2 - 9 + 8k 2 r 2 y A + A_ . 
Putting all these contributions together as dictated in (TI8~1). we end up with: 



(26) 



n+-(P) = -eVTr soft (l-x 2 ) 

9 



(8k - 8kr - 8k 2 - 5 + 8k 2 r 2 ) 2 A + A 



128k 2 



32k 3 



(k + kr- iy (r - 1) Q A_ 



(27) 



All other terms in and 11^ are treated along similar lines and we obtain for the longi- 
tudinal HTL-dressed one- loop self-energy the result: 



IT 44 (P) = -eV Tr soft (1 - x 2 ) 



9 



(r-l) 2 (l + r) Q A 



9 



32k y ' y ' ^ + 32k 

+ — (2kr -2k- l) 2 A+A' + — '— 
2k 2 K i + + 128k 2 



r + l) 2 (l-r) Q A_ 
(8k - 8kr - 8k 2 - 5 + 8A; 2 r 2 ) 2 A + a'_ 



+ (8A; 2 r 2 - 8kr - 8k 2 - 5 - 8A;) 2 A_A 

128A; 2 v 



(28) 



The imaginary part of the above expression can safely be evaluated using the Matsubara 
relation: 

ImT ^2 A e (ik , k) A v (iq , k) 

hoo r+oo 

-tt i did I dJ [1 — nf (u) — rif (u 1 )] p E (to, k) p v (— u', k) 5 (1 — u — to') . (29) 



The quantity p £ (co, k) is a spectral density and rif (u) the Fermi-Dirac distribution. Using 
this relation in (128]) and calling on ([7]) yield the result: 



7/(0) 



oo 

Arp i- /-+CJO 



e 4 T 
64tt 



dk I doo duj'5 (1 — u — u') 



e=± L 

+-k~ 2 (k + eu) 2 (k - eu) PoP ' £ 



+ ^(2u-2k-l) 2 p + p'_ 



(30) 



We have set p v = p v (cj, k) where r\ = +, — or 0, and p' v = p v (co 1 , k). Note that the integration 
over the solid-angle k is performed. 

A similar expression for the transverse damping rate 7 t (p) is longer and more tedious to 
obtain. We spare the reader the details and give the result: 



It (?) 



e 4 T 
64vr 



/p-f-OO P-f-OO 
dk doj du'[A(uj,uj',k) + C(uj,uj',k) p 2 + ...]. (31) 
J —oo J —oo 



The coefficient A (u, uj', k) is as follows: 

A = E \jU ( fc + £W ) 2 ^ - £W ) + 4 (-8w + ^8A; - 8fc 2 + 8w 2 - 5) 2 p £ 



4F 

e=± L 

+ l(2 w -2£;-l) 2 p + p'_<5, 



where 5 stands for 5 (1 — uj — uj'). The coefficient C (uj, uj', k) is: 



Pe S 



(32) 



C = E ( Cl 0eP0 + CleePe) P^ + C 1+ „p + p'_5 
e=± 

+ ^2 (CwcPo + ^WJ <9 fc p^ + C 2+ -p + d kP '_5 

e=± 

+ J] (C 30e Po + CWJ + C 3+ ^ P+ d 2 k p'_5 

e=± 

+ (CttePo + C iesPe) Ps d u S + C i+ -P+P-d u 6. 



(33) 



e=± 



In this relation, dk( w ) stands for partial d/dk(oj). The different coefficients in C above are 



8 



as follows. Those for terms without derivatives are: 

Cw£ = _ 3g(^- fc 2 ) 2 , Quj2 _ 6f(2 _ 7£k3 _ 3uj + + 3£ku 2 _ 7k 2^ 

Cue = (406ek - 758w + \2ASeku - U12k 2 + 176(k£; 3 - 1408& 4 

+ U08ek 5 - 1280A; 6 + 156w 2 + 1312w 3 + 512w 4 - 1920w 5 + 768w 6 - 864eku 2 
- 2208£; 2 w - 768eku 3 + 768ek 3 co + 384eku 4 - 29Uk 4 u + 896Jfe V + 4864Jb V 
-1792e£; 3 u; 2 - 2816fcV + 3328fcV + 143) ; 
C l+ _ = (224w -U2k- 800kco + 720k 2 + 256k 3 - 512A; 4 + 720w 2 - 1792a; 3 

+1024c/ + 3072fcw 2 - 1536£; 2 u; - 2048£;u; 3 + 1024A; 3 u; + 512JfeV - 83) . (34) 

The coefficients involving dkP £ (u/, k) are given here: 

9 2 
C 20 e = _ 4qT3 ( 2uj -^){k + euj) ( k — eu) ; 

C 2e s = —r^r flOw - 2A6k - Wekco + 88k 2 - 72u 2 + 48w 3 
480 k v 

-48A; 2 w + 7) {e8k - 8cu - 8k 2 + 8co 2 - 5) ; 

C 2+ - = — (I6fc - 16w - IQkuo + 16w 2 + 11) (2u - 2k - 1) , (35) 

15fc 

and those involving d 2 p £ (u/, k) jdk 2 write like this: 

3 2 

C 3ee = — (e8k -8cu- 8k 2 + 8co 2 - 5) 2 ; C 3+ - = ^ (2w - 2A; - l) 2 . (36) 

Finally, the coefficients involving d^S (1 — uj — co') are as follows: 

9 2 
= ^^e(ek + u) (ek-u); 

C ie£ = — (e8k -8cu- 8k 2 + 8co 2 - 5) 2 ; C 4 +- = - (2u -2k- l) 2 . (37) 
80 5 

Remember that only soft values of k, u and u' are allowed, which means we have used the 

approximation rif (co) ~ ~ (l — ^ + . . . ) . Also, the expansions in fl3]) are used to appropriate 

orders for each step of the calculations. Finally, we are still keeping oj = 1. 

Now what remains is to use the explicit expressions of the spectral densities and perform 

the integrals in ( 130]) and ( 13TI) . These are given in [l6|: 

p± (to, k) = }±{k)5 (to - u±(k)) + fc(k)6 (u + uj± (k)) + (3± (w, k) 6 (k 2 - cu 2 ) ; 
Po (co,k) = ~Q(k 2 -oo 2 ). (38) 



The residue functions i £ (k) at the fermionic poles co £ (k) and the cuts (3 £ (to, k) are known: 
k 2 - co\ (k) 



3± (*) = 
(3 ± (co,k) = 



2m) 



—m 2 (k =F to) 



2k 2 



co =f k =F 



2 2 

ffij rrif 



iTj) Qo(k/u) 



2 2 4 



(39) 



fc k 

As before, it is best to illustrate the coming steps with a generic example. Take then for 
instance the integral: 

r+oo r+oo r+oo / \ / \ 

/++ = J dk J dco J duo' f (to, k) p + (to, k) d k p + [to', kj 5 \1 — to — to' J . (40) 

Here the integral over to' is readily done with the replacement to' — 1—to. Imposing the pole 
conditions on both co and 1 — co forbids 'pole-pole' terms, so that we get: 

T - T PC 4- T CP 4- T cc - 
-<++ - 1 ++ + J ++ + I ++, 



jPC 
L ++ 



t cp 
l ++ 



T cc 
'++ 



r+oo r+oo 

/ dk dtof (oo, k) z + 5(oj - to+)d 2 k [0 + (1 - to, k) Q(k 2 - (1 - oo) 2 )] ; 

JO J-oo 

dtof (to, k) (3 + (to, k)Q(k 2 - co 2 )d 2 k [n (k) 5(1 -co- co + )] ; 

r+oo r+oo 

/ dk dcof(co,k) [3 + (co,k)0(k 2 -co 2 ) 

JO J~oc 



dk 



oo 



x^[/5 + (l- W) fc)0(fc 2 -(l- W f)]. (41) 

The usual care must be taken when manipulating the Heaviside distribution and its deriva- 
tives. For the 'pole-cut' contribution, we obtain the result: 

d k P+(w,k)\ u=1 _ u _~ 



jPC 

'++ 



J dki+f k) d 2 k (3 + (to, k) | w=1 _ w+ + 



3+/(w+,A;) 



II 



+ i 



= -0.031503. (42) 

In this expression, fc p is the solution of the equation w + (k) — 1 + fc. Also, the prime here 
means derivative with respect to fc. The contribution 'cut-pole' is given in the following 
expression: 

oo 

I c 4 = Jdk[f (co, k) P + (co, k) i" + - (2 3 ' + u/ + + u<) [f (oo, k) P + (co, k)} 



+Hu+ 2 dl[f(u,k)P + (co, fc)]] w=1 _ w+ 

d u (5 + (to, k) w=l _ u+ 



+ 



n (oo' + Yf(l-co + ,k) 



II -oo', 



+ \ 



= -2.796605. 



(43) 
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Note here that, because of the divergent behavior of the cut functions and some of their 
derivatives at k p , some of the terms, individually, are singular at this point. However, taking 
them together carefully cancels the singularities, and the final result is finite. Details have 



been presented elsewhere 



17|. Finally, the 'cut-cut' contribution writes as: 



T cc 



dk J dud k [f (u,k) (3 + (cu,k)] d k p + (l-u,k) = -0.04266. (44) 

0.5 1-fc 



It is finite, though with singularities on the light-cone handled with suitable change of 
variables and adequate algebra \x\. As we see, the integral I ++ is finite, equal to —2.870 8. 

Similar procedures are applied to all other terms with various lengths. When all are 
summed and u restored, the final results are the following: 

li (P) = ^ [1-32954 + 0(p 2 )} ; lt (p) = & — [1.32954 + 3.05041p 2 + 0(f)] . (45) 

Remember that p = p/uj®. The longitudinal and transverse one- loop HTL-dressed damping 
rates are finite, positive and equal at zero momentum. The absence of infrared sensitivity is 
related to the absence of photon propagators at one-loop HTL-dressed order, which remain 
massless beyond the electric (soft) scale. Their absence is also related to the appearance of 
the one-loop HTL-dressed photonic damping rates at the higher order e 3 (eT) instead of the 
next-to-leading order e (eT), like the fermionic damping rates for example 17| • Indeed, for 
soft ui and a/, the thermal factor (1 — rif (u) — rif (u/)) is of O (e) instead of O (1/e) if we 
had soft photons in the loops. 

It is important to stress that the results in (T4"5!) come from the calculation of the imag- 
inary part of the one-loop HTL-dressed self-energy of a soft (almost static) photon in a 
thermalized QED plasma. Loop momenta are soft. One then naturally asks: are there other 
diagrams not incorporated in the one-loop HTL summation scheme that would contribute 
to the photon damping rates at this order or smaller? A similar issue was raised previously 
in the literature, in the context of the thermal-field-theory calculation of the production 
rate of non-thermalized soft photons in a high-temperature quark-gluon plasma in thermal 
equilibrium. When calculated in the strict one-loop HTL-dressed prescription, the imagi- 
nary part of the photon self-energy Imll, to which the production rate is proportional, was 
found to be of the order e 2 g 3 T 2 , but exihibiting a colinear divergence for real light-cone 
emitted photons [h]]; g is the QCD coupling constant. This is a quite general feature of 



11 



the hard thermal loops, which become colinearly singular when some of their external soft 
momenta are put on the light cone. This problem was addressed in {20I and resolved by 
endowing the internal hard fermions with an asymptotic thermal mass moo- Single-pole col- 
inear singularities thus regularized enhance a given order by a factor In (T 2 /m 2 ^) ~ In (1/g) 
for soft moo ~ gT. However, more severely, because of double-pole colinear singularities, 
some two-loop diagrams with hard internal fermions and a soft gluon insertion that are a 
priori dismissed in the standard HTL-summation order-of-magnitude counting are found to 
contribute and dominate over the mere one- loop HTL-dressed result 21] . More precisely, 
when the double-pole colinear singularity in the two-loop diagram is regularized with moo, 
it contributes an enhancement factor T 2 /m^ D ~ 1/g 2 instead of a simple log, which means 
an order e 2 gT 2 for Imll for ligh-cone photons. This result carries through to the case of 



virtual photons, relevant to lepton pair production, but in the static 



w p h is the energy of the photon, such two- loop diagrams go to zero [211 ]. 



imit p <C w p h where 



The production rate of soft static photons has been evaluated at one-loop HTL-dressed 
order in 22] and reexamined in (23] at the two-loop level in the limit w p h <C m g <C T, where 
trig ~ gT is the gluon thermal mass. The bremsstrahlung contribution to Imll is found to 
be of the order e 2 g 2 m 2 g T/cu ph In (T 2 /m 2 g ) , of the same order as the one-loop HTL-dressed 
result in the same regime. A similar conclusion is reached for the Compton and annihilation 



24| for ^> gT, and an order e 2 g 3 at two 



processes. The same problem is reexamined in 
loops is also found. 

The photon production rate problem has been examined in different limits and approxi- 
mation schemes (25)]. In that problem, there are two independent small coupling constants, 
e and g, and three independent energy scales, T, p and u p h. Our situation is somewhat 
physically different, that of a thermalized QED plasma in which e is the only small coupling 
constant, and T and p the only two independent energy scales; the photons are thermally 
on-shell, with u p h = u)t,i (p) ~ eT. Therefore, care must be taken when drawing analogies 
between the photon production rate problem and our situation. From this work and from 



previous experience 



14j] with the case p <C w p h, we learn that colinear issues do indeed 



arise in the midst of the calculations, but never seriously enough to warrant the use of the 
asymptotic-mass improved HTL perturbation 20j so that they enhance the order of the 
contribution. It could therfore be reasonable to expect that in our situation, colinear diver- 
gences are safe enough not to enhance the order of two-loop diagrams with hard internal 
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momenta. However, it is also reasonable to expect the same two-loop diagrams to con- 
tribute to order e 5 T 2 as we cannot find a convincing order-of-magnitude counting argument 
to rule them out. But at the same time, it is also clear that only an explicit investigation 
of such diagrams will determine their order and precise contribution, as well as the nature 
of divergences they may or may not have. This task is beyond the present framework and 
is left for a future project. In light of this discussion, the results in ( 145]) must be seen as 
the contribution to the damping rates of soft photons with p <C u; p h = uj^i (p) coming from 
one-loop HTL-dressed diagrams, a contribution that is finite in the infrared and free from 
colinear divergences. 

Finally, beyond the damping rates, one would like to study the next-to-leading energies 
of the quasiparticles in HTL perturbation. For this, one needs to look at the real parts of 
the HTL-dressed next-to-leading order self-energies, something which is more complicated. 



There is already the work 



261 ] which estimated the next-to-leading gluon mass. More re- 
cently, the work |27| estimated the next-to-leading quark mass. It would be interesting to 
see if we can estimate the next-to-leading dispersion laws for slow-moving quasiparticles in 
HTL perturbation. 
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